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Abstract 



A theory of the novel spiral chaos state recently observed in Rayleigh-Benard 
convection is proposed in terms of the importance of invasive defects i.e de- 
£r} ■ fects that through their intrinsic dynamics expand to take over the system. 

o. 

The motion of the spiral defects is shown to be dominated by wave vector frus- 
tration, rather than a rotational motion driven by a vertical vorticity field. 

On. 

This leads to a continuum of spiral frequencies, and a spiral may rotate in 

o. 

either sense depending on the wave vector of its local environment. Results 

-<— > ■ 
-<— > , 

of extensive numerical work on equations modelling the convection system 

provide some confirmation of these ideas. 
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A novel spatio-temporal chaotic state in Rayleigh-Benard convection consisting of dy- 
namic spirals or targets has recently been discovered experimentally |I|, 0, and then repro- 
duced numerically first in equations that approximately model the system JJ] and then in the 
full thermally driven Navier-Stokes equations [Q. This "spiral chaos" state is particularly 
interesting because it occurs for parameter values where the familiar state of straight parallel 
convecting rolls (the "stripe" state) is known to be stable, and because the spiral geometry, 
although familiar in chemical reaction diffusion systems, was completely unexpected in this 
system where the ideal state is static. The spiral state also provides an intriguing example of 
the phenomenon of "defect chaos" where persistent dynamics is apparently associated with 
easily identified defects or coherent structures. Although numerical work has reproduced the 
experimental results with remarkable accuracy, there is still little theoretical understanding 
why the spiral chaotic state develops, or of the dynamic behavior of the spiral or target 
defects within the chaotic state. 

In this paper we introduce the idea of an invasive defect — a localized structure that 
through its intrinsic dynamics tends to grow and take over a significant portion of the 
system — as a way to understand the instability (or rather, lack of persistence) of non-ideal 
states towards the chaotic state. We analyze the dynamics of individual spirals and targets 
based on a slow- variation approach known as "phase equations" and show that "wave vector 
frustration" is crucial to understand this dynamics. Perhaps surprisingly, fluid vorticity 
plays a secondary role in the dynamics of the spirals: the spirals appear to rotate because 
of the radial motion of the rolls induced by wave vector frustration exactly as in the radial 
motion of the axisymmetric target states. This leads us to a prediction for when the spiral 
chaos state should be expected. We analyze the transition between a spiral dominated state 
and a target dominated state in terms of a novel core instability of an axisymmetric target 
to nonaxisymmetric perturbations. In addition we present results of extensive numerical 
investigations of classes of model equations that allow us to test these ideas, and for the first 
time numerically demonstrate the spiral to target transition. 

Our theoretical analysis is based on the "phase equation" || — a dynamical equation for 



a variable describing displacements of a local stripe pattern. An important feature of the 
convection system is that this dynamics is coupled to an additional long wavelength mode 
consisting of a horizontal fluid flow with nonzero mean averaged over the depth of the cell. 
Such a flow is induced by distortions of the convective rolls (e.g. curvature) and in turn 
advects the rolls. Since the mean flow must be divergence free by the near incompressibility 
of the fluid, it takes the form of a stirring motion, and may be completely described in terms 
of a vertical vorticity Q. Previous work ||, [ff] on the spiral chaos has suggested that the 
vertical vorticity resulting from the chiral nature of the spiral state may be important in the 
observed dynamics of the spirals. We find that the role of ft is more subtle than this. 

The phase variable 8(x,t) is defined as a function of time and the horizontal coordinate 
x = (x, y) so that the local wave vector is q = V#. An expansion in slow variations of a 
basic stripe pattern yields the equation |5] 

e + U-q = T- 1 {q)V-(qB(q)) (1) 

where U is the divergence free mean flow completely characterized by the vertical vorticity 
ft = z ■ V x U in turn driven by distortions of the stripes 

ft = - 1 (q)z ■ V x \qV ■ (qA 2 (q))} (2) 



In these expressions r, 7, A and B are functions of the magnitude of the local wave vector 
and depend on the control parameter as well as other fluid parameters. In (0) we have used 
the simplified form valid near threshold |J : the complicated expression valid more generally 
0, [|7| should not change the conclusions. 

From these equations we show that the spiral dynamics is induced by the radial motion 
of rolls forced by wave vector frustration rather than by a rotation of the structure driven 
by a vorticity field induced by the chiral structure. An important consequence of this is 
that a spiral of a given chiral sense may rotate in either direction. In addition the spiral 
frequencies in a system will form a continuum rather than a discrete set or unique value as 
found for chemical spirals 0]. 



This result is readily obtained by integrating ([!]) for the phase of a spiral 9 = k(r)r + 
mcf) — out over a cylindrical geometry radius R with R >> 1 for a well developed spiral: 

uj I d 2 r r(q) = / d 2 r r(q)U ■ q — / B(q)q r rd(j) (3) 

J J J r =R 

(Here r and <fi are cylindrical polar coordinates.) The second term on the right, which is also 
present in the axisymmetric target case, is zero if q = qf, the "focus selected wave number" 
|| defined by B(qj) = ||, and is nonzero only if the wave number is held away from this 
value by some other means (wave vector frustration). In this case the contribution to u is 
of order r~ 1 B'(qf)qf (q — qj)R~ x with r an average of r(q) over the cell. We can estimate 
the size of the first term by correspondingly integrating the vorticity equation |2] 

/ d 2 r [7(g)] - 1 Q = - I V • (qA 2 ) q-dl+l V • (qA 2 ) q- dl (4) 

It is readily checked that for R 3> 1 where A 2 —* const, the first term on the right hand 
side goes as R~ l . If A remained constant as r — » 0, the second term would correspondingly 
diverge as r" 1 . Instead, of course, A — > over the coherence length £ for spatial variations of 
the amplitude of the pattern. Thus the integrated vorticity is dominated by a contribution 
localized in the core region where amplitude and phase are rapidly varying, rather than 
a contribution from the slowly varying spiral arm structure at large distances. Now U ■ 
q ~ U^m/r ~ mCl^ 2 /r 2 with an average vorticity, so that the first term in (§) gives a 
contribution to the frequency u proportional to R~ 2 , smaller by the factor R~ 1 . 

The invasive nature of spirals and targets now rests on the comparison of the background 
wave number q^ in their vicinity with qf : if q^ > qf then rolls will tend to move into the 
defect center, a spiral will "unwind" and, unless there is some coherent source of these rolls 
elsewhere, the background state will tend to take over the spiral or target; on the other hand 
if qi, < qf the target or spiral will tend to invade the background. The question of what 
to assume for q^ is quite hard, and there may not be a completely general argument. We 
propose that since the background state often contains many stationary or slowly moving 
roll dislocations, an appropriate choice for qb is the value qd of the wave number at which 



an inserted dislocation has zero velocity along the rolls (climb velocity) ||, [[K]]. This leads 
to the prediction that spiral or target chaos is expected for q d < q$. 

The combined motion of a spiral and dislocation can be understood more fully by con- 
sidering the situation shown in Fig.|l| with a single dislocation terminating a spiral at radius 
R. The motion can be split into two components. There is an outward motion of the rolls, 
driven by the curvature, at a speed a(qf — q)/qR with a a positive proportionality constant 
and q the wave vector at the dislocation: this would leave the spiral as shown by the dot- 
ted line. On the other hand the dislocation will climb at a speed j3{q — q d ) moving in the 
azimuthal direction and shortening the spiral, with /3 another positive proportionality con- 
stant. A steady state solution in which the spiral rotates at constant rate and with constant 
length exists if these effects cancel i.e. if q = (aqf + f3q d )/(a + /?), a weighted mean of qf 
and q d . 

We have investigated this picture in numerical simulations of equations modelling the 
convection system as in the work of Xi et al. || . In this approach the system is described 
by an equation for a real field if)(x,t), a function of horizontal coordinates and time only, 
that represents the horizontal structure of the convection pattern. The dynamical equation 



for ip is based on the Swift-Hohenberg equation |TT| that has been much used in the study 
of stationary patterns and transient relaxation processes @ : 

^ + U ■ Vi/j = ei\) + (v 2 + l) 2 i) - gip 3 + 3(1 - g)(Vifj) 2 V 2 if) (5) 



where again U is the vorticity induced mean flow. For g = 1 fl5|) reduces to the Swift- 
Hohenberg model. The additional nonlinear term [|1^] for g =£ 1 yields a more accurate 



reproduction of the fluid stability balloon and provides an additional tuning parameter to 
investigate transitions in behavior. We have investigated various forms for the driving of 
the vorticity by distortions of the field if>: the first ("dynamic vorticity") follows Xi et al. || 

Cl - a(V 2 - c 2 )Vt = g m z- V(V 2 ?A) x Vij (6) 

where a is the fluid Prandtl number and g m determines the strength of the coupling be- 
tween the mean flow and if). We have also considered two modifications to (|6]). The first 



("passive vorticity") is to eliminate the term Q term, which together with V 2 fi is higher 
order in the gradient expansion leading to (|6|) . The second is to reduce the high wave vector 
Fourier components of the vorticity by including the "filtering operator" F 7 introduced in 



ref. [[12|].This term reduces the short wavelength components of the vorticity field, suppress- 
ing an additional instability in Eqs.(|),(H) that does not correspond to the behavior of the 
fluid system c.f. ref. [§J. 

We have fixed e = 0.7, a = 1, c 2 = 2, as in ref. |§, and have varied the vorticity coupling 
g m , the non-linearity coefficient g, and the form of the vorticity equation. Our simulations 
are done in a square geometry with periodic boundary conditions using the same numerical 



scheme as in ref. [fL3 |. A dynamic spiral or target state rapidly forms from random initial 
conditions over a wide range of parameters. We indeed see spirals rotating in either direction 
relative to their chirality, although since the unwinding sense is observed to lead to their 
destruction, the winding sense is more prevalent. 

We have studied the defect chaos state as a function of g m and g to assess the importance 
of qf and q^. Fig.|2] shows the spatial structure a time 4000 after random initial conditions 
in a size 200 x 200 using dynamic, unfiltered vorticity. Panel (a) shows the dynamic spiral 
state for g m = 50 and g = 1. As g m is reduced targets become more prevalent e.g. panel 
(b) with g m = 10. Finally as g m is reduced further, panel (c) with g m = 2, the target-spiral 
state disappears, and the pattern freezes into a disordered stripe state. If alternatively 
g is reduced changing the form of the nonlinearity, the spatial structure remains rather 
unchanged, until at low g again the spiral structures disappear, but now to a new dynamic 
state in which dislocation pairs, evident in panel (d) with g = 0.2, rapidly migrate across the 



system [14| . The quantitative aspects of the evolution of the structures with g m and g are 
shown in Fig|| Both as a function of g m and g the mean wave vector is seen to follow a trend 
consistent with a weighted mean of qf and q^, except for the lowest g m where in fact the 
spatial structure shows few spirals or targets. In addition the spiral/target state disappears 
when q<i approaches qf, although we do not have a model for which the qf and q^ lines 



actually cross. It should also be noted that the wave vector distribution is not determined 
by the stability boundaries. This is shown for example in (b) where the distribution varies 
with g m whereas the Eckhaus stability boundary does not. The irrelevance of the cross- roll 
instability is shown by similar results (not shown) for the model with filtered vorticity when 
the CR line is suppressed to larger wave vectors. 

Further insight into the spiral-target transition is given by considering the dynamics 
of a single target in an axisymmetric geometry. Fig.f| shows time independent solutions 
constrained to axisymmetry in a circular geometry with boundary conditions ip = dijj/dr = 
that act to pin the phase of the rolls at the outer boundary. As the radius is increased, the 
rolls are stretched and the value of the field at the center goes up/down, until eventually a 
new roll is nucleated at the center. For the parameters used in Fig.^ this event is hysteretic: 
if the radius is decreased the roll disappears at a smaller radius than the one at which the 
nucleation occurred. Numerical linear stability analysis shows that the portion of the curve 
with filled symbols in Fig.|] is unstable to an m = 1 symmetry breaking perturbation (the 
most unstable mode). The eigenvector of the instability, (inset to Fig.|3]), shows it to be 
localized at the core of the target, in contrast to the instability studied by Newell et al. ||TB| . 



As g m is decreased, the instability point P moves closer to the point where the hysteretic 
jump occurs, and reaches this point at g m ~ 21.5. 

In the dynamics of a target in a disordered background (e.g. Fig||) we can suppose that 
the target will approximately pass through the sequence of states represented in Fig.|]a, 
with roll radiation corresponding to increasing the system length. We cannot say exactly 
where the jump from one branch to the other, corresponding to the nucleation of a new 
roll) will occur. However we see that if the nucleation occurs sufficiently early the core 
may experience the nonaxisymmetric instability providing g m is greater than about 20. We 
might guess that this would then lead to a spiral core structure, which would further evolve 
to convert the target to a spiral. In fact our numerics of the disordered state in periodic 
boundary conditions shows that the subsequent nonlinear behavior is more complicated: in 
the target state at g m = 20 for example, we observe a complicated, almost periodic behavior 



in which the nonaxisymmetric perturbation successively grows and decays over many cycles 
(~ 15) before a new roll is finally nucleated and the structure regains the form of a target. 

In conclusion, we have proposed that the novel spiral/target chaos state may be un- 
derstood in terms of the invasive dynamics of these defects. This dynamics in which rolls 
are radiated from the center is driven by wave vector frustration. For the spirals the rota- 
tion induced by the vertical vorticity plays a secondary role in the dynamics. The invasive 
dynamics depends on the balance between the focus selected wave vector qj and the back- 
ground wave vector. The vorticity is however important in generating the spiral core, and 
the target to spiral transition may be partially understood in terms of a nonaxisymmetric 
instability of the axisymmetric core structure of a target. In addition, in the absence of the 
coupling to the vertical vorticity the focus selected wave vector gj is marginally unstable to 
a zigzag distortion, so that the vorticity may also be important in stabilizing the far field of 
the targets and spirals, as well as in the motion of the dislocations and hence the value of 
qd- Our proposal provides a criterion for when spiral chaos may be expected. The weakest 
part of our analysis is the identification of the background wave vector with q^. alternative 
possibilities might be a wave vector selected by another defect in the system that acts as a 
sink of the rolls, or a wave vector at which a roll destruction process develops, for example 
by the production of dislocation pairs or the nucleation of a mobile dislocation from the core 
of a stationary defect. 
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FIGURES 

FIG. 1. Spiral defect ending in a dislocation 

FIG. 2. Disordered states a time 4000 after random initial conditions in a system of size 
200 x 200 for dynamic unfiltered vorticity showing the variation of the patterns as g and g m 
are varied. Values are (a) g = 1, g m = 50; (b) g = 1, g m = 10; (c) g = 1, g m = 2; (d) 5 = 0.2, 
g m = 50. 

FIG. 3. Wavenumber distribution in disordered states as in fig.@ as a function of (a) g and 
(b) g m - The symbols are the mean wavevector and the error bars give the standard deviation, 
both given by a gaussian fit to the distribution. The lines show qf and the wavevectors for the 
charateristic instabilities of straight parallel rolls (E: Eckhaus; SV: skew- varicose; CR: cross-roll) 



[15], and qa given by separate numerical calculations of the velocity of climb of a pair of dislocations 



in straight roll backgrounds of various wave vectors, as in ref. flOfl . 

FIG. 4. (a) Variation of the field ip(r = 0) at the center of a target as a function of the radius 
of the system. Notice the hysteretic transition. The states represented by the filled symbols are 
unstable to an m = 1 distortion, (b) radial dependence of the eigenvector of the unstable mode at 
the point A as well as the stable eigenvector with largest eigenvalue at the corresponding point A' 
on the stable branch. Notice the localized nature of the instability. Parameters are g = 1, g m = 50, 
e = 0.7. 
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